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1  .  Introduc t  ion 

Carrier  Sense  Multiple   Access  with   Collision  Detection    (CSMA/CD)  is  a 
channel  access  protocol  for  packet  broadcasting  local  area  networks.  This 
protocol  allows  user  stations  to  contend  for  time  on  a  shared  channel  by 
multiaccessing  it  in  a  random  fashion.  In  an  attempt  to  reduce  the  frequency 
with  which  broadcast  message  packets  collide  on  the  shared  channel,  Ethernet 
(Metcalfe  and  Boggs  1976)  and  IEEE  802  (IEEE  1985)  use  a  technique  known  as  1- 
persistent  CSMA/CD  with  binary   exponential    backoff    (Goodman  et  al .  1988).  An 
alternative  technique  proposed  for  reducing  the  frequency  of  collisions  is  p- 
persistent  CSMA/CD  (see,  for  example,  Stallings  1987  and  Takagi  and  Kleinrock 
1985a);  this  protocol  is  considered  in  this  paper. 

The  dynamic  behavior  and  performance  of  unbuffered   CSMA/CD  systems  has 
been  studied  extensively  (e.g.,  Lam  1980,  Takagi  and  Kleinrock  1985a,  and 
Tobagi  1980) .  These  studies  assumed  that  there  are  an  infinite  number  of 
network  users  and  each  user  can  have  at  most  one  packet  in  the  message  queue 
at  any  time  (called  an  infinite   source  model).    This  simplifying  assumption  has 
been  relaxed  in  recent  studies  on  CSMA/CD  systems:  Apostolopoulos  and 
Protonotarios  (1986),  Goodman  et  al .  (1988),  Hastad,  Leighton  and  Rogoff 
(1987),  Park  and  Bartoszynski  (1990a),  Park  and  Bartoszynski  (1990b),  Takagi 
and  Kleinrock  (1985b),  and  Tasaka  (1986).  In  these  studies,  a  finite  number  of 
users  are  considered  to  have  message  queues  with  finite  or  infinite  capacity 
(called  a  finite   source  model    or  a  buffered   CSMA/CD  system) . 

Among  these  studies  on  buffered  CSMA/CD,  the  following  authors  considered 
p-persistent  CSMA/CD.  Under  the  assumption  of  finite  capacity  queues, 
Apostolopoulos  and  Protonotarios  (1986)  modeled  the  message  queueing  process 
as  a  two-dimensional  semi-Markov  chain.  With  the  state  space  reduced  to  a  two- 
dimensional  space,  it  was  necessary  to  employ  an  iterative  approximation 
procedure  to  obtain  the  mean  packet  delay  in  steady  state.  The  transmission 
probability  p  was  arbitrarily  set  equal  to  the  inverse  of  the  number  of  busy 
users  in  the  beginning  of  each  contention  period.  Thus  the  dependency  of  the 
mean  packet  delay  on  the  control  parameter  p  was  not  investigated.  Takagi  and 
Kleinrock  (1985b)  developed  a  stationary  joint  probability  generating  function 


(PGF  hereinafter)  of  queue  lengths  for  a  system  with  two   users,  and  obtained 
an  explicit  formula  for  the  mean  packet  delay  based  on  that  PGF.  An  exact 
joint  PGF  of  the  queue  length  vector  for  a  system  with  more  than  two  users  was 
first  developed  by  Park  and  Bartoszynski  (1990b).  Further,  they  obtained 
stability  conditions  analytically  based  on  the  PGF. 

In  implementing  p-persistent  CSMA/CD  the  major  problem  is  to  determine 
the  value  of  p  so  as  to  minimize  the  response  time  for  users  while  maintaining 
stability  of  the  system.  In  this  paper  an  approximation  procedure  is  developed 
to  compute  the  optimal  p  which  leads  to  a  stationary  operating  mode  and  at  the 
same  time  yields  the  minimum  mean  packet  delay  in  equilibrium.  An  optimization 
of  the  p-value  with  respect  to  the  response  time  has  not  been  studied  for  a  p- 
persistent  CSMA/CD  system  with  more  than  two  users,  though  it  has  been  carried 
out  for  slotted  ALOHA  by  Saadawi  and  Ephremides  (1981)  and  Sidi  and  Segall 
(1983). 

In  the  next  section,  the  p-persistent  CSMA/CD  channel  access  protocol  is 
described  and  the  queueing  dynamics  under  this  protocol  is  modeled  as  a  multi- 
dimensional semi-Markov  chain.  A  stationary  PGF  of  packet  backlog,  i.e.,  the 
sum  of  all  queue  lengths,  is  obtained  in  Section  il  based  on  the  joint  PGF  of 
the  queue  length  vector.  In  Section  4  the  Kolmogrov  forward  equation  is 
derived  applying  an  M/G/l  approximation  to  the  stationary  PGF  of  packet 
backlog,  and  an  iterative  procedure  is  developed  to  evaluate  the  mean  packet 
delay  in  equilibrium  and  optimize  p  with  respect  to  the  mean  delay.  A  numeric 
analysis  is  conducted  in  Section  5  and  the  analytic  results  are  compared  with 
simulation  results  demonstrating  the  accuracy  of  the  proposed  approximation 
scheme.  Section  6  contains  concluding  remarks. 

2.  A  Semi -Markovian  Model  of  the  Message  Queueing  Process 
A  brief  description  of  the  p-persistent  CSMA/CD  channel  access  protocol, 
which  is  sufficient  for  the  reader  to  follow  the  subsequent  analysis,  follows. 

The  system  consists  of  a  single  server  (the  channel)  and  multiple, 
infinite -capacity  queues  of  customers  (message  packets).  The  number  of  queues, 
denoted  by  m,  can  be  as  large  as  200  in  local  area  networks.  The  channel  time 


is  slotted  with  the  slot  size  being  the  maximum  propagation  delay,  and  a  slot 
size  is  chosen  as  the  unit  of  time.  Discrete  time  index  t-0,1,...  is  used  to 
denote  slot  boundaries. 

If  the  channel  is  busy  no  users  attempt  transmission  (due  to  carrier 
sensing)    until  ongoing  transmissions  are  completed  and  the  channel  becomes 
empty  again.  Once  the  channel  becomes  empty,  every  busy  user  (i.e.,  a  user 
with  a  nonempty  queue)  persistently  attempts  at  every  t   to  transmit  a  message 
packet  from  the  top  of  its  queue.  In  a  transmission  attempt  a  user  samples  a 
random  number  from  the  uniform  distribution  over  [0,1].  If  the  number  is 
smaller  than  p,    it  starts  transmission  of  a  packet;  otherwise,  it  suppresses 
the  start  of  transmission.  Thus  whenever  the  channel  is  empty  all  busy  users 
contend  to  seize  the  channel  with  an  equal  chance.  If  only  one  among  all  busy 
users  samples  a  random  number  smaller  than  p   at  time  t,    the  user  starts 
transmitting  a  packet  at  t  and  the  packet  is  successfully  transmitted.  If  two 
or  more  users  start  transmitting  packets  at  the  same  t,  packets  collide,  the 
collision  is  detected,  and  all  collided  transmissions  are  aborted  before 
transmissions  are  completed. 

It  is  assumed  that  in  each  slot,  at  most  one  packet  arrives  at  a  queue 
(which  is  reminiscent  of  Poisson  process  assumption) .  The  probability  that  a 
packet  arrives  at  a  queue  during  a  slot  is  denoted  by  X.    The  numbers  of 
packets  arriving  in  different  slots  at  the  same  user  are  independent,  and  the 
numbers  of  packet  arrivals  at  different  users  in  the  same  or  different  slots 
are  independent.  This  message  input  process  (called  a  geometric   arrival 
process)    implies  that  queueing  processes  at  different  users  are  symmetric; 
viz.,  queue  lengths  at  individual  users  are  exchangeable   random  variables 
(Billingsley  1979) . 

Following  Takagi  and  Kleinrock  (1985a) ,  the  channel  state  dynamics  are 
modeled  as  depicted  in  Figure  1  (note  that  users  3  to  m-1    in  the  figure  are 
assumed  to  be  idle  throughout  the  portion  of  time  shown  in  the  figure) .  The 
channel  state  alternates  between  idle  periods   and  busy  periods.    During  an  idle 
period  all  users  have  empty  queues,  and  during  a  busy  period  there  is  at  least 
one  busy  user.  With  the  assumed  message  input  process,  durations  of  idle 


periods  are  independently  and  geometrically  distributed;  hence  the  beginning 
of  an  idle  period  is  a  system  regenerative  point.  A  busy  period  is  divided 
into  a  number  of  sub (busy) per iods ,  each  in  turn  consisting  of  a  transmission 
delay   followed  by  a  packet    transmission  period. 


Figure  1  about  here 


During  a  transmission  delay  no  packets  are  transmitted,  even  though  the 
channel  is  empty,  since  all  busy  users  continually  sample  random  numbers 
exceeding  p.    A  transmission  delay  terminates  at  time  t   when  one  or  more  busy 
users  sample  random  numbers  smaller  than  p   and  start  packet  transmissions.  The 
duration  of  a  successful  packet  transmission  is  proportional  to  the  length  of 
the  packet  transmitted.  Packet  lengths  may  vary.  However,  for  analytic 
convenience,  we  use  the  expected  value  of  packet  lengths  and  denote  it  by  i. 
The  duration  of  an  unsuccessful  transmission,  which  is  the  time  to  detect  and 
resolve  a  collision,  is  represented  by  V  .  In  a  baseband  system  like  Ethernet, 
it  takes  at  most  two  slots  to  detect  a  collision.  Thus  2'     is  usually  much 
shorter  than  1.    A  p-persistent  CSMA/CD  system  is  completely  described  by  the 
set  of  parameters,  [m,p,  X ,  S. ,  i'  }  ,  at  least  for  the  purpose  of  our  analysis. 

Throughout  the  paper  we  use  the  following  convention  in  defining 
notations:  A  boldface  capital  letter  indicates  an  m-vector  or  a  function  with 
vector  arguments,  except  that  E  and  P  stand  for  expectation  and  probability, 
respectively.  An  italic  capital  letter  with  an  underlined  superscript,  i, 
indicates  the  ith  element  of  an  m-vector,  and  an  underlined  italic  capital 
letter  signifies  the  sum  of  all  elements  in  an  m-vector.  An  italic  capital 
letter  without  a  superscript  or  an  underscore  signifies  a  set  or  a  function 
with  scalar  arguments. 

Denoting  the  set  of  network  users  by  I={ 1 , 2 , . . , m) ,  we  define  Q  -  [Q-, 
i&I]    to  be  the  vector  of  queue  lengths  at  t,  and  Q   =  STiQt  to  be  the  packet 
backlog   at  t.    The  process  (Qr)  is  characterized  as  an  irreducible,  aperiodic 


discrete- time  Markov  chain  whose  state  space  is  the  m- dimensional  vector  space 
Q-{ [w. ,  iel}:    w.^Q,  integer). 

A  semi-Markov  chain  model  of  the  queueing  process  is  constructed  as 
follows  (which  is  more  accurate  than  the  semi -Markov  chain  model  presented  in 
Takagi  and  Kleinrock  1985b) :  The  beginning  and  the  end  of  each  idle  period  are 
regarded  as  embedded  Markov  epochs.  Within  a  busy  period,  the  boundaries  of 
subperiods  are  chosen  as  embedded  Markov  epochs .  Denoting  these  epochs  by 
k-0 , 1 , . . ,  with  k-0  at  t-0 ,  we  have  an  embedded  Markov  chain  {Q»J  which  is  a 
multidimensional  random  walk.  In  Figure  1,  the  embedded  epochs  are  shown  by  •, 
and  the  packet  backlog  at  those  epochs,  Q,,    are  illustrated. 

3.  Stationary  Probability  Generating  Function  of  Packet  Backlog 

The  major  difficulty  in  analyzing  the  queueing  process  (Qk)    lies  in  that 
queues  at  different  users  interfere  with  each  other  through  the  shared 
channel.  This  problem  of  interacting  multiple   queues    cannot  be  solved  using 
classical  queueing  theory.  There  are  no  analytic  results  in  the  literature  for 
evaluating  the  mean  packet  delay  in  such  a  multidimensional  queueing  system, 
and  even  approximations  are  difficult  to  obtain  when  m>   2    (see  Kleinrock  and 
Yemini  1980) .  In  this  paper  we  base  the  derivation  of  the  mean  packet  delay  in 
equilibrium  on  the  stationary  PGF  of  {<2tJ  .  By  the  aggregation  of  (QjJ  into  a 
one-dimensional  process  (Q/J  we,  of  course,  lose  information  about  behaviors 
of  individual  queues.  The  approximation  needed  to  overcome  this  loss  of 
information  can,  however,  be  localized  in  the  delay  analysis  (as  will  be  shown 
in  Section  4)  so  as  to  generate  accurate  results  (verified  by  simulation  in 
Section  5) . 

To  describe  the  process  (2;.}  we  first  define  the  following  random 
variables.  Notations  for  each  element  in  a  random  vector  and  the  sum  of  all 
elements  in  a  random  vector  are  defined  according  to  the  convention  mentioned 
earlier. 
X(n)   :  numbers  of  packet  arrivals  at  respective  queues  over  an  n-slot  period. 


Z,     :  numbers  of  packets  that  are  successfully  transmitted  from  respective 

queues  at  the  end  of  the  subperiod  which  begins  at  the  kth  epoch;  Note 
that  either  Z,=0  (0-vector)  or  Zr  -1  for  some  i  and  Zt   -0  for  all  ;>i . 

R,  :  duration  (in  number  of  slots)  of  the  transmission  delay  in  the  sub- 

period  beginning  at  the  kth  epoch. 

L,  :  duration  (in  number  of  slots)  of  the  packet  transmission  period  in  the 


Uk 
Bk 


Bk<V 


subperiod  beginning  at  the  kth  epoch. 

■  { iGl :  (hr   >1 }  ,  i.e.,  the  set  of  busy  users  at  the  kth  epoch. 

■  {iel:    Qg-0 ,X-(R. )>1 ) ,  i.e.,  the  set  of  users  who  are  idle  at  the  kth 
epoch,  but  become  busy  by  the  end  of  the  transmission  delay  in  the 
subperiod  beginning  at  the  kth  epoch. 

The  process  iQi.)    can  now  be  recursively  defined  as 

2k+l"  Qk^k     where  ^k=  £(1)  if  2k=0;  ^k*  ^W'^k  othervise-    (1) 


In  this  equation  £(1)  represents  random  variable  £(1)  with  the  constraint 
£(1)>  1.  Under  the  geometric  input  process,  X(l)    has  a  binomial  distribution 
with  parameters  m   and  A.  Therefore,  if  0.^=0    (i.e.  ,  if  the  kth  epoch  is  the 
beginning  of  an  idle  period)  ,  fir., -i  follows  a  truncated  binomial  distribution 
given  by 

?lQk+l=q\Qk=0)    =  \^\qa-l)m'q/(l-(l-\)m),      1<  q<  m.  (2) 

Note  that  2r.+i  given  S^'O  is  the  packet  backlog  at  the  start  of  a  busy  period. 
Next,  if  Q,>    1,  the  random  drift  U.     is  the  sum  of  random  variables  X(R,), 
X(L,)    and  -Z,  .    Unfortunately,  the  latter  variables  are  stochastically 
dependent  on  each  other  in  a  complex  way.  First,  L,     (and  hence  X(L.))    depends 
on  Z,.    Namely,  L,=2'    if  Z,=0   and  L.— i   otherwise  (i.e.,  if  Z, -1) .  Random 
variable  Z,     in  turn  depends  on  |B.|  and  \B!(R,)\    since  the  success  of  a  packet 
transmission  depends  on  the  number  of  busy  users  at  the  end  of  the  transmis- 
sion delay,  which  is  the  sum  of  \B,\    and  \B',(R,)\.    The  latter,  \B',(R,)\, 
depends  on  X(R,).    Further,  R.    and  X(R,)    are  jointly  dependent  on  \3,   | ,  and 


finally  \B,\    is  uniquely  determined  given  Q^,  but  stochastically  dependent  on 
Q,.    These  stochastic  dependencies  are  summarized  in  Lemma  1  below.  A  detailed 
analysis  of  these  stochastic  dependencies  and  the  resultant  conditional  PGF  of 
Q,  1  given  Q^-l  can  be  found  in  Park  and  Bartoszynski  (1990b) .  A  derivation  of 
the  latter  PGF  and  the  proof  of  Lemma  1  are  summarized  in  Appendix  I  at  the 
end  of  this  paper. 


Lemma  1 :  For  any  k  with  Qk>0 . 
z(U.k\Qk)-  tt\l(Rk)+ZZ 


ULk) 


Zk 


-Ik 


\B'k(Rk)\ 


Rk,l(Rk) 


<2,  • 


An  optimization  of  the  p-value  with  respect  to  the  mean  packet  delay  has 
been  studied  for  slotted  ALOHA,  as  mentioned  earlier  in  Section  1.  Compared 
with  slotted  ALOHA,  the  delay  analysis  and  optimization  of  the  p-value  is 
significantly  more  complex  for  p-persistent  CSMA/CD.  The  source  of  additional 
complexity  is  the  stochastic  dependencies  among  the  random  variables  as 
described  in  Lemma  1.  These  stochastic  dependencies  stem  from  the  carrier 
sensing   and  collision   detection   mechanisms  which  have  been  added  to  the 
slotted  ALOHA  protocol  in  devising  p-persistent  CSMA/CD  (see  Szpankowski  1988 
and  Tsybakov  and  Mikhailov  1980  for  a  Markov  chain  model  of  slotted  ALOHA  with 
buffered  users) .  In  this  paper  the  stochastic  dependencies  among  L, ,  Z,   ,  R,  , 
K(Rk)    and  | B,   |  are  exactly   incorporated  in  the  delay  analysis.  An 
approximation  is  introduced  to  simplify  the  dependency  of  \B,   |  on  0. . 

The  process  iQ.k)    is  a  random  walk  with  random  drifts  iHj.)  •  The  process  is 
uniformly   downward   bounded   meaning  that  P(L/,<-l)-0  for  all  k.    It  also 
possesses  an  important  property  called  the  bounded  homogeneity ,  which  is 
described  in  the  following  Lemma  (see  Appendix  I  for  the  proof): 

Lemma  2 :  For  any  k,    k'  ,    Qk  and  Qk,    such    that    \B,\  =  \B,,  \  ,    P(^-u|2^)- 
?(Uk,=u\Qk,). 


Define  F,(s)    to  be  the  PGF  of  Q, .  In  view  of  (1)  we  write 


Fk+i(s)-  Fk+i^l2k-0)p(V0)  +  ^k+i(s\ak>o)?(Qk>o) 


(3) 


Letting  H(s):~   F^.+1  (s  1 2^.-0 )  .  we  obtain  from  Equation  (2) 


H(s)-    {(Xs+l-X)m-(l-X)m)/{l-(l-\)m) . 


(4) 


By  Lemmata  1  and  2  we  may  write 


Fk+1(s|<2k>0)-  E 


s  kEE<U  k 


|SJ 


2^>0 


,      °o    min(q,ra)  <=° 


P(0,>0)  s  ,s         ^  , 

ViK   '  q=l     0=1   u=-l 


^  su  P(t7k=u|  |Sk|=b;  P(|Sk|-b|2k-q)  P(2k-q).   (5) 


Define  Gk(s  \  \B.  \  ~b)  :  -  z"=_  ,sUP(17,=u  |  |  S.  |  =5)  ,  i.e.,  the  conditional  PGF  of 

U,    given  |3,|«b.  By  Lemma  2,  for  any  b^O ,  G(s|b):=  G,(s||S,|=b)  is  independent 
of  k   and  Q,  .  By  setting  s—  =s  for  all  iEl  in  Equation  (A4)  in  Appendix  I,  we 
get 

\s 


G(s\b)> 


m-b    r 

m-b 
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I  J  J 


i(l-p)(As+l-A)-(l-A)j 
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"1(As  +  l-A)mi(b+j)p(l-p)5+^"1  +(As+l-A)'"r  (l-(l-p)5+-'-(b+j)p(l-p)fa+J"1 


V7-0 


(J)(-i)J"n 


l-(l-p)    (1-A)      (As+l-A) 


(6) 


Note  that  index  j    in  Equation  (6)  represents  the  number  \B',(Rk)\. 

By  substituting  (6)  into  (5),  (4) -(5)  into  (3),  and  then  by  passing  k  to 
=o  in  (3),  the  stationary  PGF  of  { Q,  )  ,  denoted  by  F(s)  ,  is  written  as 


<*>    min(  q,m) 

F(s)~  H(s)?(QrO)   +  S  sq         2    G(s|b)P(|S|-b|2-g)P(2-q).      (7) 

q=l     b-1 


As  in  (7)  we  simply  drop  the  time  index  k   from  an  argument  to  indicate  its 
quantity  in  steady  state. 

4.  Delay  Analysis 

In  the  sequel,  we  write  F  -*    (c)  to  denote  the  Jth  partial  derivative  of 
F(s)    with  respect  to  s,  evaluated  at  s=c .    We  need  to  evaluate  F(   (1)  to 
obtain  the  mean  packet  backlog  at  an  embedded  epoch  in  steady  state,  which  in 
turn  is  used  to  evaluate  the  mean  packet  delay.  Analytic  evaluations  of  F(s) 
and  its  derivatives  are  difficult  to  obtain  mainly  because  the  conditional 
probability  P(|B|=b|o=q)  cannot  be  evaluated  explicitly  in  the  presence  of 
interactions  among  queues.  Therefore,  we  introduce  an  M/G/l  approximation  in 
formulating  this  probability. 

In  this  approximation,  we  assume  that  in  steady  state  the  message  queue 
at  a  user  behaves  like  a  discrete  M/G/l  queue  independently   of  other  queues. 
The  queue  at  a  user  in  the  CSMA/CD  system  indeed  has  a  single  server  (the 
channel)  and  geometric  arrivals  (a  discrete  version  of  Poisson  arrivals)  of 
message  packets.  But  how  about  the  service  time?  The  service  time  for  a  packet 
begins  at  the  point  when  the  packet  is  promoted  to  the  top  of  its  queue.  The 
service  time  is  defined,  in  the  context  of  an  M/G/l  system,  to  consist  of:  (i) 
a  number  of  sub(busy)periods ,  since  that  beginning  point,  during  each  of  which 
the  packet  is  not  transmitted  (due  to  sampling  a  random  number  greater  than  p) 
or  is  involved  in  a  collision;  plus  (ii)  the  subperiod  in  which  the  packet  is 
finally  transmitted  with  success.  This  service  time  is  actually  affected  by 
other  busy  users'  contention  for  the  server.  Furthermore,  the  set  of  busy 
users  may  change  during  the  service  time  of  a  particular  packet.  The 
assumption  of  independence  (or  "no  interference")  among  the  queues  in  steady 
state  is  thus  introduced  so  that  service  times  are  i.i.d.  in  steady  state. 

With  this  assumption,  we  now  may  write 


P(|B|-b)-  fyeba-e)m-b 


(8) 


where  8    is  defined  to  be  the  steady-state  probability  of  a  user  having  a  non- 
empty queue  at  an  embedded  epoch. 

The  no   interference   assumption  for  the  steady-state  behavior  of  multiple 
queues  should  lead  to  a  closer  approximation  when  it  is  applied  to  p- 
persistent  CSMA/CD  than  when  applied  to  slotted  ALOHA.  This  is  because  the 
former  suffers  less  from  interference  (viz.,  collisions)  than  the  latter  due 
to  additional  controls,  carrier  sensing   and  collision   detection .    The 
approximation  in  (8)  is  motivated  by  the  results  in  Sidi  and  Segall  (1983) 
where  an  M/M/l  approximation  is  applied  to  the  slotted  ALOHA  demonstrating 
good  accuracy  over  a  wide  range  of  system  parameters.  It  should  be  noted  that 
the  approximation  is  limited  only  to  computing  8 ,  and  the  rest  of  delay 
analysis  exactly  reflects  the  dynamic  behavior  of  interacting  multiple  queues. 

Incorporating  (8)  into  (7)  and  using  the  equation  F(0)-  P(2-0),  we  obtain 


F(s\d)~  H(s)F(0\d)    +  V(s\9){F(s\8)-F(0\8)) , 


(9) 


m 


;here 


V(s\8): 


G(s\b) 


5=1 


(?) 


8ba-9)m-b 


(10) 


Taking  the  first  derivatives  of  both  sides  of  (9)  yields 


F(1)(s\d)~      H(1) (s)F(0\8)    +  V(1) (s\d){F(s\d)-F(0\d) )     /(l-V(s\d)).       (11) 


/ 


By  evaluating  Equation  (11)  at  s=l  and  using  the  condition  F(1|0)-  l/(l|#)-  1, 
the  equilibrium  differential  equation  (i.e.,  Kolmogrov  forward  equation  in 
equilibrium:  see  Kleinrock  1975)  can  be  derived  as 


0  =  H(1)  (l)F(0\8)    +  \/<rI  ^  ( 1 1  ^  )  {  1  -F(0  |  ^  )  }  , 


(12) 
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which  can  be  rewritten  as 

F(O|0)-  V(1) (l\e)/{V(1) (l\0)-H(1) (l)) .  (13) 

Here  JT  '(1)  and  V(  '(l|0)  are  the  expected  drift  of  packet  backlog  during  an 
idle  period  and  that  during  a  sub(busy)period,  respectively.  From  (4)  we  get 
H(1)(l)    -  mX/{l-(l-\)m) .    V(1)(l\d)    is  obtained  simply  by  replacing  G(s\b)    by 
G^\l\b)    in  the  RHS  of  Equation  (10).  G^  \l\b)    is  the  expected  drift  of 
packet  backlog  over  a  subperiod  beginning  with  b   busy  users,  and  is  derived  in 
Appendix  II  using  Equation  (6) . 

In  light  of  (12),  the  condition  v'   ^(l|0)<  0  must  be  satisfied  for  the 
system  to  achieve  equilibrium.  In  Equation  (10)  we  see  that  V(    '(l\d)    is  the 
marginal  expectation  of  G(  '(l|b)  over  b=l,..,m,  which  in  turn  is  equivalent 
to  E(£/,  \Q.t>0)  .    Thus  the  equilibrium  equation  (12)  implies  that  a  necessary 
condition  for  stability  of  p-persistent  CSMA/CD  (i.e.,  ergodicicy   of  (Qr))  is 
approximately   given  by  E(l\  |(2,>0)<  0.  Furthermore,  by  Pakes'  (1969)  lemma,  the 
latter  condition  is  also  sufficient  for  ergodicit;-.-  of  the  one -dimensional 
random  walk  {Qj.)  ■    It  can  be  easily  proved  that  ( Q_ )  is  ergodic  if  and  only  if 
{Q,  }    is  ergodic.  Therefore,  the  condition  E(t/,|(2.  >0)<  0  can  be  regarded  as  a 
sufficient  and  necessary  condition  for  stability.  However,  a  theoretical  proof 
of  this  assertion  is  difficult  and  has  not  been  established  in  the  literature. 

The  reader  is  referred  to  Park  and  Bartoszynski  (1990b)  for  a  theoretical 
analysis  of  stability  conditions  for  p-persistent  CSMA/CD.  Some  of  the  main 
results  in  that  paper  are  summarized  in  the  following  for  later  reference  in 
Section  5.  The  conjecture  below  has  been  numerically  established. 

Lemma  3:  (Qc)  is   recurrent   if  G(1)  (1 1  b)-E(t/fc|  \B,  \  =b)  <  0  for  all    b-l,..,m,    and 
only   if   mA[^+(l-p)m/{l-(l-p)n7}]  <  1. 

Conjecture :  Given   system  parameters    {m,\ , 2 , £' } ,    for  any  b    there   exist   PL(b) 
and   pu(b)  such  that   G(l)(l\b)<   0  if  and   only   if  PL(b)<  p<   Py(b) • 
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Further,     (Qc)  is   recurrent    if   PLd)^  P^  PyCxn)  since  pL(b)<  p    (b ' )  and 
p,,(b)<  p,,(&')  for  any  b>  5 ' . 

Recall  that  9    is  the  probability  of  a  user  being  busy  in  steady  state, 
and  that  F(0\d)    is  the  conditional  probability  of  all  users  being  idle  in 
steady  state  given  the  value  of  9  .    It  follows  that,  under  the  assumption  of 
independent  queues,  we  must  have  1-9  —  F(0\d)    '     .    We  thus  compute  the  value  of 
9    iteratively  using  Equation  (13)  until  the  following  convergence  criterion  is 
met : 

\9    -    (l-F(0\9)1/m) I  <  €.  (14) 


The  function  f{9)\-  9+F(0\$)  '  -1  is  monotone  increasing  in  9  and  changes  sign 
once  in  the  range  (0,1)  of  9.  Therefore,  a  bisecting  search  method  can  quickly 
find  the  value  of  9    satisfying  (14)  even  for  a  extremely  small  value  of  t. 

Let  9  be  the  value  of  9  satisfying  (14).  We  can  now  evaluate  F  (l\B)  by 
applying  L'Hospital's  rule  to  (11)  and  then  replacing  the  RHS  of  Equation  (13) 
for  F(0\9) : 


.(1) 


(i|0>- 


■  CD 


(2) 


(1)  7u;(l|9)  -  H{ZJ\  I)    V{1J  (l\9) 


(2) 


»d) 


/ 


/ 


2    V(1)(l\9) 


(1)(l\d)    -    H(1)(l) 


(15) 


(2) 


(2) 


In  this  equation,  tfuy(l)=  m(ra-l) Az/{ l-( l-X)m)    from  (4),  and  V{ZJ(l\9)    is 
obtained  simply  by  replacing  G(s\b)    by  C        (l\b)    in  the  RHS  of  Equation  (10). 
A  derivation  of  G(   (l|b)  is  provided  in  Appendix  II. 

Recall  that  F        (l\9)    is  the  mean  packet  backlog  observed  only  at  the 
embedded  Markov  epoch.  In  order  to  obtain  the  steady- state  mean  packet  delay 
observed  at  an  arbitrary  slot  boundary  t,  we  need  to  evaluate  the  mean  packet 
backlog  at  an  arbitrary  t    in  steady  state  (see  Takagi  and  Kleinrock  1985b  for 
a  similar  treatment).  If  t  is  in  an  idle  period,  the  backlog  is  zero.  The 
probability  chat  an  arbitrary  t    is  found  in  a  busy  period  when  the  system  is 
in  steady  state  is  given  by  { 1  - F( 0 | 9 ) } .  Letting  Q   denote  the  mean  packet 
backlog  at  an  arbitrary  t  in  steady  state,  we  have 
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2-  (i-F(o|$)}s/r 


(16) 


where  T   is  the  expected  duration  of  a  sub(busy)period  in  steady  state,  and  S 
is  the  expected  backlog  accumulation  over  a  subperiod  in  steady  state. 

First,  to  evaluate  T,    we  note  that  the  duration  of  a  subperiod  is  i.i.d. 
for  all  subperiods  starting  with  the  same  number  of  busy  users,  due  to  the 
bounded  homogeneity    (see  Lemma  2).  Letting  6(b)  denote  the  expected  duration 
of  a  subperiod  starting  with  b  busy  users,  we  may  write  T-   E?, 6 (b)P( |B |-b) . 
Furthermore,  due  to  regeneracy   of  the  queueing  process  the  duration  of  a 
subperiod  starting  with  b    (>0)  busy  users  converges  a.s.  to  E(i?,  +L,  |  \B,  |-b) 
provided  that  the  process  is  ergodic.  Consequently,  for  any  q>0  we  have 


6(b)=   E(i?1+L1|  |B1|=b)=  EJi^+E 


h 


Rl,\B[(Rl)\ 


li'l-bj. 


(17) 


reflecting  the  stochastic  dependencies  among  R,   ,  L,  and  Z,     (see  Lemma  1). 

Define  T(s|b)  to  be  the  conditional  PGF  of  (R^+L-,)    given  |S,|-5>0.  This 
PGF  and  5(b)-  r^  '(l|b)  are  derived  based  on  (17)  in  Appendix  III.  Using  (A9) 
in  the  Appendix  and  Equation  (8)  we  get 


m 


b=l 


r(1)a\b)[™} 


lb(l_d)m-b 


(18) 


Turning  to  S    in  (16),  let  us  first  define  /3(q)  as  the  expected  backlog 
accumulation  over  a  subperiod  beginning  with  q   outstanding  packets  so  that  S— 
S  ^P(q)?(Q-q)  ■    To  derive  /3(q)  for  q>0 ,  we  define  £(y\q,b)    to  be  the  total 
time  (in  number  of  slots)  that  l£L)  spends  in  state  y   during  the  subperiod 
which  starts  with  q  packets  waiting  at  b    (>1)  busy  users.  Due  to  bounded 
homogeneity ,    £(y\q,b)    is  i.i.d.  for  all  subperiods  that  start  with  q  waiting 
packets  and  b   busy  users.  Further,  it  converges  a.s.  to  its  expected  value  if 
the  queueing  process  is  ergodic.  Therefore,  for  any  q>0  we  have 
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«  min(q,/n) 

/Kg)-  ^     2    y  E£(y|q,b)  P(|B|-b|2-q).  (19) 

y-q    b-1 


The  process  (Cj  }  is  non-decreasing  during  a  subperiod  and  the  packet  arrival 
process  is  time  homogeneous.  Consequently,  the  average  backlog  at  an  arbitrary 
slot  boundary  in  the  subperiod  which  starts  with  q   waiting  packets  and  b  busy 
users,  or  E™  yE£  (y  |  q  ,  b)/<5  (b)  ,  is  estimated  by  q+(  .  5)mXS  (b)  .  We  thus  can  write 

<d  <=°  min(q,m) 

S-      ^(q)P(2-q)~   )      )     6(b){q+(.5)mX6(b))    P( | S | -b |£-q)P(£-q) .    (20) 
q-1  q-1    b-1 

From  (9)  we  see  that  2°°  ,  qP(2=q)=  d[  V(s  I  0  )  {F(s  I  5  )  -F(0  I  0  )  }  ]/ds  \       ,-  F^(l|0)  + 

q-1  |  s-1 

If'  '(l|0){l-F(O|*)}.  Incorporating  the  latter  and  (8)  into  (20)  we  get 


m 
S-  t\f(1) (l\$)+V(1) (l\$){l-F(0\ti))      +  (.5)mA  \  T(1)  (1 1  b)  2  H  5b(l-  d)m'b  .   (21) 

b-1 

Finally,  using  Little  (1961) 's  theorem  the  mean  packet  delay  in  steady 
state,  D,     is  given  by 

D-  Q/m\-    (1-F(0| 9))S/mXT,  (22) 

The  mean  packet  delay,  D,  includes  (i)  the  elapsed  time  from  the  arrival  of  a 
packet  until  the  packet  is  promoted  to  the  top  of  the  queue,  (ii)  a  number  of 
sub(busy)periods  during  each  of  which  the  packet  is  not  transmitted  or  is 
involved  in  a  collision,  and  (iii)  the  subperiod  in  which  the  packet  is 
successfully  transmitted. 


5.  Numeric  and  Simulation  Experiments 
The  algorithm  to  compute  D   given  parameters  (p ,m,X , £ , 2' }    was  encoded  in 

Fortran  77;  it  first  finds  d    and  F(0|#)  satisfying  the  condition  (14)  with 
e-l.E-10,  and  then  compute  D   using  Equations  (15)  through  (22). 
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Computational  results  are  shown  in  Tables  1  and  2,  respectively,  for  two 
different  cases.  The  two  cases  have  the  same  values  for  m    (-50)  and  i.'    (-3), 
but  different  values  for  A  and  2.    The  latter  parameters  were  set  to  1.6E-4  and 
75,  respectively,  in  Case  1  to  represent  a  high  utilization  of  the  channel 
capacity;  and  to  1.2E-4  and  25,  respectively,  in  Case  2  to  represent  a  low 
utilization.  Note  that  by  Lemma  3  an  upper  bound  on  the  channel   capacity 
(i.e. ,  the  maximum  achievable  throughput  rate  in  equilibrium)  can  be  given  by 
1/i .  The  ratio  of  the  total  input  rate,  mX ,  to  1/i  is  62%  in  Case  1  and  15%  in 
Case  2. 


Tables  1  and  2  about  here 


In  view  of  the  conjecture  in  Section  4,  the  range  [p, (1) ,p„(m) ]    of  p   that 
guarantees  stability  is  provided  in  each  table.  Among  the  steady-state 
measures  reported  in  the  tables,  the  new  notation,  B    ,  signifies  the  set  (b: 
G^   ' (l\b)>   0).  The  p-values,  at  which  v'   '(l|$)<  0  and  B+  is  nonempty,  yield  a 
finite  mean  packet  delay,  though  not  satisfying  the  sufficient  condition  for 
stability  given  in  Lemma  3.  It  should  be  noted,  however,  that  the  delay  at 
p=p.y(/n)  is  close  to  the  minimum  delay:  in  Case  1  the  minimum  delay  was  128.6 
at  the  optimal  p-value  .16,  while  the  delay  at  p=p,,( 50)=-.  08  was  140.4, 
exhibiting  a  9%  gap  from  the  minimum;  in  Case  2  the  gap  was  only  2%.  It  is 
also  notable  that  the  expected  number  of  busy  users  in  steady  state,  or  E(|B|) 
(=  md ) ,  was  1.5  and  1.1,  respectively,  in  Cases  1  and  2  when  p  was  set  at  the 
optimum.  This  indicates  that  there  is  little  interference  among  queues  in 
steady  state  if  p  is  optimized  with  respect  to  the  mean  packet  delay,  which 
justifies  the  M/G/l  approximation.  Indeed,  the  ratio  of  collisions  to  all 
packet  transmissions  observed  in  simulation  runs  was  about  5%  and  1%  in  Cases 
1  and  2,  respectively,  with  the  optimal  p-values. 

The  simulation  model  was  written  in  SLAM  II  (Pritsker  1984)  and  run  on  an 
IBM  3033AP  mainframe.  In  the  simulation  model,  we  generated  packet  arrivals 
for  each  user  so  that  the  interarrival  times  follow  an  exponential 
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distribution  with  intensity  A.  It  is  well  known  that  the  geometric  distribu- 
tion of  interarrival  times  can  be  closely  approximated  by  the  exponential 
distribution  when  A  is  small  (Feller  1968) .  This  exponential  approximation 
tremendously  reduced  the  computing  time.  Approximately  5,000  packet  transmis- 
sions were  collected  to  estimate  the  mean  packet  delay  for  each  p-value.  The 
first  500  observations  were  discarded  to  eliminate  the  initial  transient 
period  (see  Wilson  and  Pritsker  1978,  Schruben  and  Goldsman  1985). 

The  simulation  results  are  shown  in  Figure  2  along  with  the  numerical 
approximation  results.  The  solid  line  and  dotted  line  show  approximation 
results  for  Cases  1  and  2,  respectively.  Simulation  results  for  Cases  1  and  2 
are  plotted  using  symbols  '  +  '  and  '*',  respectively,  and  show  good  agreement 
with  approximation  results.  The  delay  (versus  p)    curve  is  U-shaped  and  the 
minimum  delay  is  attained  in  the  neighborhood  of  the  maximum  p-value  that 
leads  to  equilibrium  (viz.,  near  the  bottom  right  corner  of  the  U  curve).  A 
similar,  yet  different,  pattern  of  the  delay  curve  was  found  for  slotted  ALOHA 
systems  by  Sidi  and  Segall  (1983).  The  CSMA/CD  achieves  close-minimal  mean 
packet  delays  over  a  much  wider  range  of  the  p-value  than  the  slotted  ALOHA, 
presumably  due  to  the  additional  medium  access  control  mechanisms  —  carrier 
sensing  and  collision  detection. 


Figure  2  about  here 


6  .  Cone lus  ion 
A  multidimensional  queueing  process  with  interactions  among  individual 
queues  arises  in  many  computer  and  communication  systems  such  as  coupled 
processors  (Szpankowski  1988) ,  ALOHA  satellite  communication  (Sidi  and  Segall 
1983)  and  CSMA/CD  local  area  networks.  However,  there  are  no  analytic  results 
for  computing  the  mean  queueing  delay  in  such  a  system  with  more  than  two 
queues . 
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In  this  paper  we  have  shown  that  the  mean  packet  delay  in  CSMA/CD 
networks  can  be  obtained  quite  accurately  using  an  M/G/l  approximation.  The 
success  of  this  approximation  method  is  ascribed  to  the  following:  (i)  The 
message  queueing  process  has  the  bounded  homogeneity  property  (Lemma  2).  With 
this  property,  after  introducing  the  approximate  distribution  of  the  number  of 
busy  users  (Equation  (8)),  we  could  simplify  the  stationary  PGF  of  the  packet 
backlog  (Equation  (7))  into  a  product  form  (Equation  (9)).  (Note  that  in 
Equation  (9),  Z*  ,  P(fi-q)  in  Equation  (7)  is  factored  out  as  (F(s | d) -F(0 | 8) ) . ) 
(ii)  The  delay  analysis  incorporates  an  exact  PGF  for  the  random  drift  of  the 
packet  backlog  over  a  period  between  two  consecutive  embedded  Markov  epochs 
(Equations  (2),  (6)  and  Lemma  1).  (iii)  Finally,  the  interference  among  queues 
is  actually  minimal  in  steady  state  if  the  control  parameter  p  is  optimized. 

A  practical  and  easy  way  to  determine  the  p-value  would  be  to  set  it 
equal  to  p,,(m)  ;    this  guarantees  stability  (by  Lemma  3  and  Conjecture)  and  at 
the  same  time  yields  a  close-minimal  mean  packet  delay  in  equilibrium  (as 
numerically  demonstrated  in  Section  5).  According  to  Conjecture,  for  given 
system  parameters  {m,X ,  £ ,  2'  )  ,    p,j(/n)  is  obtained  by  solving  <j   ^(l|m)=  0  for  p. 

Substituting  m   for  b    in  Equation  (A5)  in  Appendix  II,  we  obtain 


G(1)(l\n)- 


m 


1-d-p) 


m 


\a-p)m+{m\2-\)p{l-p)m-l^X2'  [l-a-p)m-mPa-p)m'1) 


■    (23) 


Rewriting  (23)  we  see  that  G(  '(l|m)«  0  is  equivalent  to 


/(q)  :-  -cxq   +c7q        +c3=-  0, 


(24) 


where  q-l-p,  c:-  c2+A(i'-l),  c2-  m\{2- 2' } -I   and  c3=  \2'  .    The  equation  /(q)  has 

■k 

the  global  minimum  at  q   ={m-\)c2/vac-i  .    It  has  exactly  two  solutions  (q^.q^) 
such  that  0<  qx<   q  <   q2<   1  provided  that  m\2<   1  (which  is  a  necessary 
condition  for  stability  by  Lemma  3)  and  f(q    )<  0.  A  bisecting  search  over  the 
range  (0,q  )  of  q   can  quickly  find  ql  ,    and  then  we  have   p,.(m)-  1-Qi- 
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Our  approach  --  viz.,  incorporating  an  exact  transient  analysis  of 
multiple  interacting  queues  into  an  M/G/l  approximation  for  evaluating  average 
delay  in  steady  state  --  can  be  employed  to  analyze  other  CSMA/CD  protocols. 
It  would  be  interesting  to  attempt  applying  this  approach  to  the  analysis  of 
new  variants  of  CSMA/CD  that  are  designed  for  high-speed  broadcast  bus 
networks  (e.g.,  p . -persistent   protocol   by  Mukherjee  and  Meditch  (1988),  LCSHA- 
CD   by  Maxemchuk  (1988),  and  a  modified  CSMA/CD  by  Lin  and  Sousa  (1990)). 


REFERENCES 
Apostolopoulos ,  T.  K. ,  and  Protonotarios ,  E.  N.  1986.  Queueing  Analysis  of 

Buffered  CSMA/CD  Protocols.  IEEE  Trans.    Commun.    COM- 34,  898-905. 
Billingsley,  P.  1979.  Probability  and  Measure.    Wiley,  New  York. 
Feller,  W.  1968.  An   Introduction    to  Probability  Theory  and  Its  Applications , 

Vol.    I.    third  ed. ,  Wiley,  New  York. 
Goodman,  J.,  Greenberg,  A.  G.,  Madras,  N.,  and  March,  P.  1988.  Stability  of 

Binary  Exponential  Backoff.  J.    Assoc.    Comput .    Mach .    35,  579-602. 
Hastad,  J.,  Leighton,  T. ,  and  Rogoff,  B.  1987.  Analysis  of  Backoff  Protocols 

for  Multiple  Access  Channels.  Proc.    ACM  Symposium  on   Theory  of  Computing , 

241-253. 
IEEE.  1985.  Carrier  Sense  Multiple  Access  with  Collision   Detection    (CSMA/CD) 

Access  Method  and  Physical   Layer  Specifications .    ANSI/IEEE  Std  802.3. 
Kleinrock,  L.  1975.  Queueing  Systems,    Vol.    I:    Theory.    Wiley,  New  York. 
Kleinrock,  L. ,  and  Yemini ,  Y.  1980.  Interfering  Queueing  Processes  in  Packet- 
Switched  Broadcast  Communication.  Proc.    IFIP   Congr . ,  557-562. 
Lam,  S.  S.  1980.  A  Carrier  Sense  Multiple  Access  Protocol  for  Local  Networks. 

Comput.    Networks   4,  21-32. 
Lin,  B.,  and  Sousa,  E.  S.  1990.  A  Modified  CSMA-CD  Protocol  for  High  Speed 

Channels.  Proc.    IEEE   INFOCOM ,  253-259. 
Little,  J.  D.  C.  1961.  A  Proof  for  the  Queueing  Formula  L=XW .    Opns .    Res.    9, 

383-387. 
Maxemchuk,  N.  F.  1988.  Twelve  Random  Access  Strategies  for  Fiber  Optic 

Networks.  IEEE  Trans.    Commun.    COM- 36,  942-9  50. 
Metcalfe,  R.  M. ,  and  Boggs ,  D.  R.  1976.  Ethernet:  Distributed  Packet  Switching 

for  Local  Computer  Networks.  Commun.    Assoc.    Comput.    Mach.    19,  395-404. 
Mukherjee,  B.,  and  Meditch,  J.  S.  1988.  The  p . -Persistent  Protocol  for 

Unidirectional  Broadcast  Bus  Network.  IEEE  Trans.    Commun.    COM- 36,  1277- 

1286. 
Pakes ,  A.  G.  1969.  Some  Conditions  for  Ergodicity  and  Recurrence  of  Markov 

Chains.  Opns.    Res.    17,  1058-1061. 


19 


Park,  J.  S.,  and  Bartoszynski ,  R.  1990a.  Queueing  Analysis  for  Buffered 
Persistent  CSMA  Local  Networks.  Infor   28,  190-200. 

9 

Park,  J.  S.,  and  Bartoszynski,  R.  1990b.  Stability  of  p-Persistent  CSMA/CD. 

Working  Paper  No. 90-15,  College  of  Business  Administration,  The  University 

of  Iowa,  (also  In  First  ORSA  Telecommunications   SIG  Conference ,  143-149, 

Boca  Raton,  Florida,  1990;  Submitted  to  J.    Assoc.    Comput .    Mach.) 
Pritsker,  A.  A.  B.  1984.  Introduction    to  Simulation  and  SLAM  II.    second  ed. , 

Halsted  Press,  New  York. 
Saadawi,  T.  W. ,  and  Ephremides ,  A.  1981.  Analysis,  Stability  and  Optimization 

of  Slotted  ALOHA  with  a  Finite  Number  of  Buffered  Users.  IEEE  Trans. 

Automat.    Contr.    26,  680-689. 
Schruben,  L. ,  and  Goldsraan,  D.  1985.  Initialization  Effects  in  Computer 

Simulation  Experiments.  Technical  Report  #594,  School  of  OR  &  IE,  Cornell 

Univ. ,  Ithaca,  NY. 
Sidi,  M. ,  and  Segall,  A.  1983.  Two  Interfering  Queues  in  Packet-Radio 

Networks.  IEEE  Trans .    Commun .    C0M-31,  123-129. 
Szpankowski,  W.  1988.  Stability  Conditions  for  Multidimensional  Queueing 

Systems  with  Computer  Applications.  Opns .    Res.    36,  944-957. 
Stallings,  W.  1987.  Local   Networks .    MacMillan,  New  York. 
Takagi,  H.,  and  Kleinrock,  L.  1985a.  Throughput  Analysis  for  Persistent  CSMA 

Systems.  IEEE  Trans.    Commun.    C0M-33,  627-638. 
Takagi,  H.,  and  Kleinrock,  L.  1985b.  Mean  Packet  Queueing  Delay  in  a  Buffered 

Two-User  CSMA/CD  System.  IEEE  Trans .    Commun.    C0M-33,  1136-1139. 
Tasaka,  S.  1986.  Dynamic  Behavior  of  a  CSMA-CD  System  with  a  Finite  Population 

of  Buffered  Users.  IEEE  Trans.    Commun.    C0M-34,  576-586. 
Tobagi,  F.  A.,  and  Hunt,  V.  B.  1980.  Performance  Analysis  of  Carrier  Sense 

Multiple  Access  with  Collision  Detection.  Comput.    Networks   4,  245-259. 
Tsybakov,  B.  S.  and  Mikhailov,  V.  A.  1980.  Ergodicity  of  a  Slotted  ALOHA 

System.  Prob.    Infor.    Transmission   15,  301-312  (translated  from  Russian 

original  in  Prob.    Peredaci    Infor.    15,  72-89.) 
Wilson,  J.  R.,  and  Pritsker,  A.  A.  B.  1978.  A  Survey  of  Research  on  Simulation 

Startup  Problem.  Simulation    31,  55-58. 

20 


Appendix  I.  The  Conditional  Joint  PGF  of  Q,.,-,    Given  0,*0 
The  joint  probability  of  R^-r   and  X(r)-[xi , iel:    0<x^<r]  conditional  on  Q, 
is  given  by 

r  iB/ci  +  i5k(r)hr r(r^xi        r-xn 

P(i?k-r,X(r)-[xi,ieI]|Qk)-  |l-(l-p)   *    *    J|  |  |  (1-p)  r[  J  J  X    l(l-X)        *j   X 


r 

{  TT   a-*)r}  ITT   5  U-A)w""1Mi-p)r"w 


r-w 


x.-l     r-v-x.+K 
A     (1-A)      2   J.  (Al) 


In  this  equation,  w   represents  the  slot  in  which  the  first  packet  since  the 
beginning  of  the  kth  subperiod  arrives  at  user  i  (who  had  no  waiting  packets 
in  the  beginning  of  the  subperiod) .  It  should  be  noted  that  B!{r)    is  uniquely 
determined  given  B,  and  X(r) ,  and  B,  is  uniquely  determined  given  Q,  .  Equation 
(Al)  thus  shows  that  the  joint  distribution  of  R,    and  X(R, )  depends  only  on 
B,   ,  apart  from  the  system  parameters  p,  X,    and  I.  That  is,  the  joint 
distribution  is  identical  for  any  k   and  k'    with  B,=B.   ,  ,  even  if  Q.*Q.  ,. 

The  probability  of  having  a  collision  in  the  subperiod  beginning  at  the 
/cth  epoch,  given  Rk=r ,    is 

r  |B,|  +  |B'(r)|  |B,  |  +  |B'(r)|-K   , 

P(Zk=0|Bk,B^(r))=  |l-(l-p)  K  K  -(\Bk\  +  \B'k(r)\)pa-p)      K  K  j/ 


{(1-d-p) 


\Bk\  +  \B'k(r)\ 


} 


(A2) 


The  probability  of  having  a  successful  transmission  in  the  subperiod  beginning 
at  the  fcth  epoch,  conditional  on  Rn=r,  is  simply  given  by  P(Z.  5*0  |b,  ,B',  (r)  )=-  1- 
P(Zk=0|B^,B^(r))  .  As  such,  given  R,=r,  Z,  depends  only  on  B,  and  BUr)  .  Recall 
that  B',(R,)    is  uniquely  determined  by  B,  and  X(£,  )  and  that  R,    and  X(i?.  )  are 

jointly  dependent  on  B,  .  Therefore,  Z,  is  transitively  dependent  only  on  B,  . 
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Since  X(L. )  depends  only  on  Z^(  it  follows  that  the  random  drift  U,  depends  on 
Q,  only  through  B.  .    We  thus  have  proved  both  Lemmata  1  and  2. 

Let  Ffc+1  ( s  |  Q^°)  .  where  s-fs^.iel],  be  the  conditional  PGF  of  Q,    given 
Qk*0,  and  G(s|S):-  G^(s|S^-S)  be  the  conditional  PGF  of  U.  given  B ,-B  which  is 
independent  of  k   and  Q,  (by  Lemma  2) .  Then  we  have 


Fk+1(s|Qk*0)=  E 


f             Q±    (             Uj- 

1 

" 

[iel           [iel 

-kJ 

Qk-0 

4 

1    m 

Y  \  G(s\B)\Yk(s:sr0  m$B)-      ^  F^Cs  :  s_- ,  -0  ,  *_.-<)  Vi^B)  j  .  (A3) 

k*0)b=l  VBd:  IB  1=5  i  'eB 


P(Q^O),^ 


Using  (Al)  and  (A2),  function  G(s|B)  is  evaluated  as 


G(s|B)=  E<     I  I  s.  E 


iel 


s  .    E 

i 


X%fc) 


•-iti 


s  . 

i 


Z,. 


^.X(^) 


Vs}- 


co    r 


r=0  x.=0  Viel 

i 


iel 


J  z.=0  Viel 


X-(Lk) 


iel 


Z.=[z.,iel' 


-z  . 

1 


iel 


X  ?\zk=[zL,ieI)    Rk-r,X(r)-[xi>iEl],Bk-BJ  P^fc-rfX(r)-[xi,iEl]  Bk-sj 


m-  b 


j-0  VB'Cl-B: |B' |=j 


}      {   II   As 


(l-p)(Asi+l-A)-(l-A) 


l> 


X  «[[n(s)p(i-P)b+J'1.    ^     s:1]  +  $(s){i-(i-P)5+J'-(b+j)p(i-p)b+-'-1}j.  x 


ieBuB  ' 


(-DJ'n 
.•2-0   VAcB  '  .  \A\=n 


b+n  /-,    s  .m-b-n 


l-(l-p)"Tii(l-A)"'  u   ll       ||  (As  +1-A) 

ieBuA 


-1 


(A4) 
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£' 


where  II(s)-    (As.+l-A)   and  *(s)-    (As.+l-A)   are  the  PGF's  of  the 
i€l         l  iel         l 

duration  of  a  successful  and  an  unsuccessful  transmission  period, 
respectively.  A  detailed  derivation  of  Equation  (A4)  can  be  found  in  Park  and 
Bartoszynski  (1990b). 


Appendix  II.  Derivatives  of  the  PGF  C(s\b)    Evaluated  at  s-1 
Taking  the  first  and  second  derivatives  of  G(s|b),  which  is  given  in 
Equation  (6),  with  respect  to  s  and  evaluating  the  derivatives  at  s-1  we 
obtain 


c(1)a\b)- 


m-b 


j-0 


m-b 

'  x  , 

J 

0 

[  J  J 

U-pJ 

hi 

l-(l-p)b+J'}jp(l-A)/(p-A)  +    (Ami-l)a  + 


+  Ami' (l-(l-p)b+^-a) 


f  Xfill- a-p)b+j) 


and 


(A5) 


c(2)a\b)=  ^ 


m-b 

1 

J-0 


m-b 

'    X    ' 

J 

r 

[  J   J 

U-pJ 

L 

a{A2mi(mi-l)-2Ami+2)  +   A2mi  '  (mi  '  - 1)  { 1-  (1  -p)b+i '-a  ) 


+  2 


(Ami-l)a  +  Ami' { 1 - (1 -p) b+J -a 


X<f>   +   jp(l-A)^/(p-A) 


+ 


+  {l-(l-p)b+j) 


AZ(2*+0')  +  2jpA(l-A)0/(p-A) 


+ 


+  j[l-a-p)b+J)4> 


j-1  +  2jA(l-p)/(p-A)  +  (j+l)A2(l-p)2/(p-A)2 


,  (A6) 


where  a:-  (6+j)p(l-p) 


b+j-l 


;  p:-   l-(l-p)    (1-A)      ;  yj>: 


n=0 


pn)(-D>^; 


4>:=     2  [i](-T-)j'n(b+n)(l-p)/p2;    *':-  S    £>]  ( -  l)j"n(b+n)  (b+n-1)  (l-/3)//32  ; 
r.=0  n=0 
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and  n : - 


2  (^)(-i)j"n(b+n)2(1-^)2/^3 

n-0 


Appendix  III.  Derivation  of  the  PGF  T(s|b)  and  Its  First  Derivative 
In  this  Appendix,  a,    0   and  if>   defined  in  Appendix  II  are  used  again.  Based 
on  Equation  (17)  we  can  write 


T(s|b)-  z\s    l   E 


h 


R1,\B{(R1)\ 


i.l-bj- 


°°   m-b  r      ,\    1 

m-b 


r=0   j-0 


I  J  J 


z-0 


-l=z 


p[zi=z||B|(r)|=j,|S1|-b]p[i?k-r,|B^(r)|-j||B1|-b]. 


oo     m-b  r 

m-b 


r-0    j-0 


I  J  J 


£      £ ' . ,     , ,   N b+j   , 
s  a   +  s   {1- (1-p)  J  -a) 


a-x)w-lxa-P)r-w 


Lw-l 


(1-A) 


r(zn-b-j)(1_p)rb  x 


(A7) 


Evaluating  the  summation  over  w  and  then  the  infinite  series  over  r,  we  obtain 


m-b 


m-  b 


r(s|b)-  ) 

j=o^ 


U-pJ 


s*a  +  si'tl-Cl-p)b+J-a}]{  >  [^](-l)J'"n/{l-s(l-^))].,(A8) 
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ra;(1 
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>  j   . 


L^-pJ 


£q  +  V  (l-(l-p)5+J-a 


+  {l-(l-P)b+J}  X 
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Table  1.  Steady  State  Measures  for  Case  1  where  m-50 , i-75 , 1 '-3 , A-l . 6E-4 

(pL(l)-1.6344E-2,  pu(50)-8.4378E-2) 

p  B+       E(|B|)  V(1)(l\ti)      F(1)(l\9)  T  S  D 


1.144E-2 

1 

* 

co 

1.308E-2 

1 

2.2886 

-.1067 

3.9976 

98.32 

427.39 

491.15 

1.634E-2 

- 

2.0754 

-.1370 

2.9448 

91.44 

296.99 

357.26 

1.961E-2 

- 

1.9466 

-.1598 

2.4314 

86.31 

233.08 

291.19 

2.288E-2 

- 

1.8611 

-.1773 

2.1347 

82.41 

195.91 

252.57 

2.768E-2 

- 

1.7778 

-.1964 

1.8762 

78.13 

163.28 

218.49 

3.902E-2 

- 

1.6721 

-.2242 

1.5871 

71.82 

126.27 

179.65 

5.036E-2 

- 

1.6194 

-.2398 

1.4565 

68.11 

109.17 

161.73 

6.170E-2 

- 

1.5887 

- .2494 

1.3832 

65.63 

99.31 

151.51 

7.594E-2 

- 

1.5655 

-.2570 

1.3273 

63.43 

91.55 

143.64 

8.438E-2 

- 

1.5563 

- .2601 

1.3047 

62.42 

88.29 

140.43 

1.097E-1 

39 

- 

50 

1.5405 

- .2655 

1.2611 

60.14 

81.63 

134.19 

1.434E-1 

29 

- 

50 

1.5337 

-.2679 

1.2312 

58.00 

76.36 

129.89 

1.603E-1 

26 

- 

50 

1.5336 

-.2679 

1.2223 

57.14 

74.49 

128.62 

1.727E-1 

25 

- 

50 

* 

co 

*   The  system  is  unstable  with  v'   ' (l\d)>   0. 
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Table  2.  Steady  State  Measures  for  Case  2  where  ro-50 , £-25 , i '-3 , A-l . 2E-4 

(pL(l)-4.8303E-2,  pu(50)-1.0722E-l) 


p 

B+ 

E(|B|) 

V(1)(l\&) 

F(1)(l\$) 

T 

S 

D 

2.415E-3 

1,2 

* 

ao 

2.898E-3 

1 

2.4953 

- .0841 

6.5053 

139.23 

964.57 

1065.35 

3.381E-3 

1 

2.1872 

-.1200 

4.1352 

128.33 

577.65 

670.06 

4.347E-3 

1 

1.8287 

- .1843 

2.3321 

109.55 

284.79 

366.04 

4.830E-3 

- 

1.7190 

-  .2111 

1.9372 

101.95 

220.65 

297.99 

7.245E-3 

- 

1.4373 

- .3040 

1.1881 

76.62 

97.77 

163.21 

2.190E-2 

- 

1.1502 

- .4556 

.7534 

37.61 

22.98 

70.03 

3.896E-2 

- 

1.0991 

-  .4920 

.7062 

28.19 

14.25 

56.52 

5.361E-2 

- 

1.0831 

-.5042 

.6928 

24.75 

11.66 

52.25 

7.309E-2 

- 

1.0733 

- .5119 

.6846 

22.24 

9.97 

49.45 

1.072E-1 

- 

1.0677 

-.5163 

.6791 

19.96 

8.59 

47.32 

1.608E-1 

33  - 

50 

1.0691 

-  .5152 

.6775 

18.21 

7.66 

46.33 

1.823E-1 

29  - 

50 

1.0712 

-.5135 

.6778 

17.77 

7.46 

46.27 

2.037E-1 

26  - 

50 

1.0738 

-  .5115 

.6784 

17.39 

7.30 

46.28 

2.144E-1 

24  - 

50 

* 

00 

*   The  system  is  unstable  with  V^   '(l|0)>  0. 
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Figure  2.  Comparison  between  Approximation  and  Simulation:  D   versus  p 
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